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Much  work  has  been  done  on  screening  procedures  under  the 
assumption  of  a  bivariate  normal  distribution.  However,  very 
little  effort  has  been  expended  on  data  which  are  from  a  trun¬ 
cated  bivariate  normal  distribution.  Methods  are  developed  for 
a  screening  procedure  to  increase  acceptable  product  from  a 
truncated  distribution.  An  acceptance  criterion  on  a  linear 
combination  of  the  largest  order  statistics  from  a  truncated 
normal  population  with  a  given  truncation  point  is  discussed. 

This  paper  also  uses  the  approximate  distribution  of  the  sample 
correlation  coefficient  in  random  samples  of  any  size  drawn  from 
a  singly  truncated  bivariate  normal  distribution  to  obtain  a 
lower  confidence  limit  on  the  population  correlation  coefficient 
p .  The  screening  procedure  discussed  here  is  based  on  knowledge 
of  the  truncation  point,  the  sample  size  and  the  lower  confidence 
limit  for  p . 

KEY  WORDS  AND  PHRASES:  Singly  truncated  bivariate  normal  distri¬ 
bution;  performance  variable;  screening 
variable;  acceptance  sampling;  truncated 


normal  distribution. 


1 .  INTRODUCTION 


In  developing  screening  procedures,  many  methods  can  be  utilized, 
depending  on  the  data  we  have  and  the  nature  of  the  problem.  Most  of 
the  previous  work  done  in  the  area  of  screening  procedures,  e.g,,  Owen- 
Boddie(1976) ,  Owen  et.  al.(1975  and  1977)  ,  is  based  on  data  from  a 
bivariate  normal  distribution,  which  is  utilized  to  calculate  the  pro¬ 
portion  successful  after  selection.  In  this  paper  we  will  consider 
screening  procedures  when  the  data  available  are  from  a  truncated 
bivariate  normal  distribution. 

A  performance  variable  with  a  one-sided  specification  cannot  be 
measured  directly,  but  a  related  variable  (called  a  screening  variable) 
can  be  measured.  In  the  language  of  acceptance  sampling  the  performance 
variable  may  be  lifetime  or  some  other  variable  for  which  the  act  of 
measuring  would  degrade  the  item.  If  the  quality  control  engineer  only 
keeps  records  of  those  values  of  the  performance  variable  and  of  the 
screening  variable  for  the  acceptable  product,  then  the  values  of  the 
performance  variable  must  exceed  some  lower  limit,  say  w^.  Hence,  past 
data  may  only  be  available  on  a  singly  truncated  bivariate  normal  distri¬ 
bution. 

Let  Y  be  a  future  performance  variable  and  X  be  a  future  screening 

variable  having  a  joint  bivariate  normal  distribution  with  parameters 

(u  ,a2,u  ,a2;p).  Let  W  be  the  past  screening  variable  which  exceeds 
x  x  y  y 

wn,  i.e.,  W  has  a  truncated  normal  distribution  with  parameter  ( u  ,cr2). 

Let  Z  be  the  past  performance  variable  and  then  (W,Z)  follow  a  singly 
truncated  bivariate  normal  distribution  where  only  W  is  truncated.  A 
one-sided  lower  specification  limit,  L,  is  given  on  the  future  performance 
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variable  Y,  i.e.,  all  items  with  Y  values  above  L  are  acceptable  and 

those  with  Y  values  below  L  are  not  acceptable.  Suppose  that  the  pro~ 

portion  of  acceptable  items  in  the  future  is  y  before  screening.  That 

is,  the  proportion  above  the  lower  specification  before  screening  is  y. 

The  screening  procedure  is  set  up  to  raise  the  proportion  successful 

from  y  to  5,  where  6  >  y.  Our  procedure  will  be  to  accept  all  items 

for  which  X  >.  aW,  N+  bw  ,  where  b  =  l-a  and  W,  .  is  the  largest  order 
(n)  o  (n) 

statistic  from  a  truncated  normal  distribution  (px,cr^)  based  on  a  sample 

of  size  n.  After  screening  we  want  to  be  100n %  sure  that 

Pviv  IT  {Y  %  L| X  i  aW,  .  +  bw  }  £  d,  i.e.,  the  proportion  of  Y's 

Y|X*  (n),  P  W  ° 

greater  than  L  is  at  least  6  in  the  screened  population.  The  reason 

for  using  aW^  +  bwQ  instead  of  a  linear  combination  of  the  sample 
n 

mean  W  =  £  W./n  and  sample  standard  deviation  S  is  that  after  trun- 

n  i=l  1  W 

eating  a  normal  distribution,  many  nice  properties  no  longer  exist  and 

to  our  knowledge  no  manageable  expression  for  the  joint  density  function 

of  W  and  STT  has  yet  been  derived, 
n  W 


2 .  PROBABILITY  EXPRESSION 

Consider  the  case  where  the  parameters  ^x>ax>^y  an< 3  °y  are 
known  and  p  is  unknown.  Let  P{Y  5  L)  =  y  and  P{X  >,  wq)  *  p  be  given. 


We  make  the  transformations  -  (X  -  y^)/^,  =  (Y  -  y  )/a^, 

U  ■  (W,  -  w  )/o  ,  V  *  Z,  -  aU.  Let  (L  -  y  )/<J  =  -K  and 

(n)  o  x’  1  y  y  y 

(wQ  -  Vx)/ox  *  -K  .  The  problem  then  becomes  one  of  finding  "a"  such  that 
P[P{Z^  >  -K^jv  >  ~Kp}  i  <S]“  n  where  the  outer  probability  is  with 
respect  to  the  estimator  of  p  and  the  inner  probability  is  the  conditional 
normal  given  the  screening  procedure  and  the  correlation.  Then 


J 


P{Z„  >  -K  |v  >  -K  }  can  be  written  as 

2  y 1  p 


af  [G(x-K  )-l+p]nG'  ( -ax  +  K  )  G 

A  =  —fl _ E _ _ _ £ 


fapx  +  K  y  _  pKp 


l  ^1- 


dx 


nP[G(x-K  )-l  +  p]n-1G'(x  -  K  )  G  (-ax  +  K  )dx 
op  p  p 


(2.1) 


i  2  * 

where  G’(z)  *  (2ir)  T  exp(-z  /2)  and  G(z)  »  f  G’(t)dt  for-®  <  z  <  ®, 

*JX> 

are  the  univariate  normal  density  and  cumulative  distribution,  respectively. 

Theorem  2.1  Under  the  assumptions  given  in  sections  1  and  2, 

A  =  P_,iv  rt  (Y  >,  l|x  5  aW.  +  bw  }  is  an  increasing  function  of  p. 

Y|X,  x,  p  (n)  o 

The  proof  of  this  theorem  is  obtained  by  showing  th^c  the  numerator 

of  equation  (2.1)  is  a  monotonically  increasing  function  of  p  since  the 

denominator  is  free  of  p  and  positive. 

Suppose  that  we  want  to  be  100n%  sure  that 

Py|x  W  p{Y  5  k  x  2  aW(n)  +  bwo}  '  tbat  is»  that  >,  6}  *  n. 

This  is  equivalent  to  P{p  >,  p*}  >,  n  for  some  p*.  Our  goal  is  to  find 

a  and  b(=  -a  +  1)  such  that 

PV|V  U  *  LIX  '  +  bw  ^  ■  lS* 

Y ' X»"(n) (n)  o 

In  order  to  solve  for  a,  p*  has  to  be  known.  Once  p*  is  known,  the 


problem  which  remains  is  to  solve  the  following  equation  for  a: 

dx 


af  [G(x  -  K  )  +  p  -  l]nG'(-ax  +  K  )G 
OP  p 


ap*x  +  K  -  p*K  I 
Y  P 


/I  -  p*z 


nf  [G(x  -  K  )  +  p  -  l]n-1G'(x  -  K  )G(-ax  +  K  )dx 
OP  P  P 


6  . 


(2.2) 


Solutions  to  this  equation  will  be  discussed  in  section  5. 
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3.  DISTRIBUTION  OF  THE  SAMPLE  CORRELATION  COEFFICIENT 


Let  (W^jZ^Xi  -  1,  . . . ,  n)  be  a  random  sample  of  size  n  from  the 


past  record  which  follows  a  singly  truncated  bivariate  normal  distribution. 


where  ^  wq,  i  *  1,  n.  In  the  following  discussion  we  consider 


the  standardized  singly  truncated  bivariate  normal  distribution,  in  which 
only  W  is  truncated.  Extension  to  the  non-standardized  case  is  straight¬ 


forward  by  using  new  variables  (u  +  a  W,  u  +  a  Z)  instead  of  (W,Z) .  Let 
7  e  x  x  ’  y  y 


Rt  be  the  sample  correlation  coefficient.  Let  F^C.jp)  and  FR(.;p)  be  the 
respective  distribution  functions  of  R^,  and  R.  Applying  equation  (32)  of 
Gayen  (1951),  the  c.d.f.  of  R^  is  given  by 


V/'V  i  FR(r;oT>  +  (C1L41  +  C2L61> 


+  (ClL42  +  C2L62)7rFR(r;PT)+C2L63frVr=‘>T)'  0'» 

3pT  3pt 


where  C 


C1  - 

8n(n  + 1)  and 

II 

H 

hJ 

3pT(X40  +  "W 

II 

<M 

PT(X40  +  X04)  - 

L61  = 

-15pT(X32Q  +  Xq23, 

+  18(X30X21  +  X( 

L623 

-94(X30  +  X03> 

+  30pt(X3qX21  + 

L63“ 

“PT(X30  +  X03) 

+  ^  V21  + 

n  -2 


2  12n(n+l) (n+3) ’ 

-  4(X31  +  X13)  +  2ptX22, 

4pT(X3i  +  X13)  +  2(2  + 

)  -  9Pt^X21  +  X12)  +  UX12X21 

03X12)> 


2)X 

r  22* 


3t)(X21  +  X12^  +  36pTX21X12 


03  12 

•  3pT(4 

i  X  )  . 
03  12''  ’ 


2Ul  2 
3T  X21 


& 


r  2i  12 


p^,  is  the  correlation  coefficient  between  W  and  Z  and  the  X 


are  the  semi-invariants  of  the  singly  truncated  bivariate  normal 


A 


distribution. 


The  results  of  the  theoretical  distribution  of  R^,  were  checked  by 
comparison  with  a  Monte  Carlo  simulation.  For  each  of  the  sample  sizes 
n  *  3,  15,  50,  truncation  points  w^  *  0.  (-.5)  -3.  and  p*-,90(,10)  .90, 

4000  values  of  R^,  were  generated.  We  made  comparisons  between  the 
empirical  and  the  theoretical  c.d.f.  of  R^,  based  upon  the  Kolmogorov- 
Smirnov  test  and  conclude  that  the  approximation  holds  well. 

4.  A  LOWER  CONFIDENCE  LIMIT  ON  p 
As  we  have  seen  from  Section  2,  a  100ri%  lower  confidence  limit  p* 
on  p  is  required  for  our  acceptance  criterion  X  aW(n)  +  (1  -  a)wQ 

Let  F  (.;p)  be  the  c.d.f.  of  R^  By  the  probability  integral  transfor- 

X 

mation  theorem,  F^CR^p)  follows  a  uniform  (0,1)  distribution.  It 

follows  that  P{Fd 
*T 

P{g(p)  S  0  }  =  n.  From  the  inequality  g(n)  $  0,  we  would  like  to  get  an 
inequality  p  >  p*  so  that  p*  is  a  100n%  lower  confidence  limit  on  p .  To 
do  this,  we  need  R^,  n  and  w^.  Since  the  exact  distribution  of  R^,  is  not 
known,  we  will  use  the  approximate  distribution  of  R^  given  by  equation 
(3.1).  It  can  be  shown  that  the  function  g(p)  is  a  decreasing  function 
of  p.  For  each  given  confidence  coefficient  n,  p*  can  be  obtained  by 
examining  the  root  of  the  equation  g(p)  =  0. 

Since  -1  <  p*  <  1,  we  can  use  the  IMSL  (1979)  subroutine  ZFALSE,  i.e., 
the  false  position  method,  to  find  the  root  of  g(p)  ■  0.  Tables  1-4  give 
the  result  of  this  computation  for  n  ■  15,  50  and  R^,  =  .40,  .50. 


(R^p)  $  n>  -  n.  Let  g(p)  *  FR^(RT;p)  -  n;  then 
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5 .  SCREENING  CRITERION 

For  given  n,  n,  Y>  <$(  >  y)  and  WQ(  =  -Kp),  a  100n%  lower  confidence 

JL 

limit  p  can  be  found  using  the  method  in  Section  4.  Thus  a  is  the  only 

unknown  in  equation  (2.2).  We  use  the  Gauss  Laguerre  quadrature  to  approxi- 

r"  -x 

mate  integrals  of  the  form  JQe  h(x)  dx  with  h(x)  a  function  of  x.  The 

solutions  can  be  found  iteratively  by  using  an  algorithm  due  to  Miller 

(1956),  which  finds  the  zeros  of  nonlinear  functions.  The  screening 

criterion  is  X  5  aW,  .+  (1  -  a)w  . 

(n)  o 


6.  EXAMPLE 

Let  Y  be  the  performance  variable  for  some  device  which  is  expensive 
to  measure.  Assume  X  and  Y  have  a  joint  standard  bivariate  normal 
distribution  with  unknown  correlation  coefficient.  A  sample  of  size  n  =  50 
is  taken  from  the  past  record  in  which  all  the  performance  scores  are  at 
least  0,  and  the  highest  performance  score  in  this  sample  is  W^  ~  2.6. 

The  sample  correlation  coefficient  R^,  is  found  to  be  .50.  Suppose  we  want 
to  be  95%  sure  that  the  proportion  of  acceptable  items  will  be  raised  from 
Y  =  .70  to  6  *  .90  after  screening.  We  use  linear  interpolation  in  Table  4 
and  find  p*  =>  .45814.  Then  we  compute  a  from  equation  (2.2)  and  it  is 
a  =  .37280.  Thus  our  screening  criterion  is  to  accept  all  items  for  which 
X  >  .96928. 


7.  CONCLUSION 


It  has  always  been  necessary  to  solve  screening  problems  first  by 
assuming  all  parameters  known.  Then  estimates  of  the  parameters  based  on 
a  training  set  are  used.  In  this  paper  we  have  assumed  that  y^,  o^,  y^,  a 
and  w^,  the  truncation  point  were  known.  Obviously  there  still  exists  the 
unsolved  problem  of  what  to  do  when  any  of  these  parameters  are  unknown. 
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TABLE  1 

CONFIDENCE  COEFFICIENT  ETA  WHEN  N  =  15  AND  RT=  .40 


W0  0.000  -0.500  -1.000  -1.500  -2.000  -2.500  -3.000 
RHO* 

.65  .99930  .99973  .99990  .99996  .99998  .99998  .99999 
.60  .99871  .99942  .99975  .99989  .99993  .99995  .9-9995 
.55  .99783  .99890  .99947  .99972  .99983  .99986  .99987 
.50  .99658  .99807  .99896  .99940  .99960  .99967  .99969 
.45  .99488  .99685  .99812  .99882  .99916  .99929  .99933 

.40  .99262  .99509  .99682  .99785  .99837  .99859  .99866 
.35  .98971  .99266  .99487  .99628  .99705  .99739  .99751 
.30  .98600  .98937  .99205  .99388  .99493  .99542  .99560 
.25  .98136  .98502  .98809  .99031  .99166  .99231  .99255 
.20  .97562  .97937  .98266  .98517  .98676  .98757  .98788 

.15  .96859  .97213  .97538  .97796  .97967  .98057  .98092 
.10  .96006  .96298  .96577  .96809  .96967  .97053  .97088 
.05  .94977  .95156  .95333  .95485  .95592  .95651  .95676 
.00  .93745  .93745  .93745  .93745  .93745  .93745  .93745 
.05  .92276  .92017  .91747  .91503  .91322  .91218  .91173 

.10  .90533  .89921  .89269  .88667  .88214  .87950  .87835 
.15  .88474  .87398  .86234  .85145  .84316  .83829  .83616 
.20  .86049  .84387  .82568  .80852  .79539  .78765  .78426 
.25  .83205  .80823  .78200  .75720  .73823  .72709  .72221 
.30  .79880  .76639  .73071  .69710  .67157  .65669  .65023 

.35  .76008  .71775  .67144  .62828  .59596  .57738  .56941 
.40  .71518  .66179  .60415  .55147  .51285  .49109  .48194 
.45  .66340  .59021  .52937  .46821  .42472  .40089  .39109 
.50  .60409  .52709  .44833  .38108  .33517  .31089  .30121 
.55  .53681  .44910  .36327  .29376  .24874  .22598  .21725 

.60  .46151  .36579  .27756  .21090  .17053  .15122  .14415 
.65  .37891  .28000  .19582  .13769  .10533  .09087  .08588 
.70  .29103  .19616  .12357  .07892  .05655  .04734  .04439 
.75  .20200  .12045  .06632  .03764  .02501  .02030  .01892 
.80  .11897  .06007  .02780  .01367  .00834  .00657  .00611 

.85  .05228  .02099  .00776  .00320  .00177  .00135  .00125 
.90  .01252  .00362  .00098  .00031  .00015  .00011  .00010 
.95  .00052  .00008  .00001  .00000  .00000  .00000  .00000 
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TABLE  2 

CONFIDENCE  COEFFICIENT  ETA  WHEN  N=  15  AND  RT= 


WO  0.000  -0.500  -1.000  -1.500  -2.000  -2.500  -3 
RHO* 


-.60 

.99968 

.99987  .99995  .99998 

.99999 

.99999 

-.55 

.99943 

.99973  .99988  .99995 

.99997 

.99997 

-.50 

.99906 

.99951  .99975  .99987 

.99992 

.99993 

-.45 

.99853 

.99915  .99953  .99972 

.99981 

.99984 

-.40 

.99780 

.99861  .99915  .99945 

.99960 

.99966 

-.35 

.99682 

.99783  .99855  .99899 

.99922 

.99932 

-.30 

.99553 

.99672  .99763  .99823 

.99857 

.99872 

-.25 

.99384 

.99519  .99628  .99705 

.99751 

.99772 

-.20 

.99167 

.99311  .99434  .99525 

.99582 

.99610 

-.15 

.98891 

.99032  .99159  .99259 

.99323 

.99357 

-.10 

.98543 

.98664  .98779  .98872 

.98936 

.98970 

-.05 

.98105 

.98182  .98259  .98324 

.98369 

.98394 

.00 

.97558 

.97558  .97558  .97558 

.97558 

.97558 

.05 

.96877 

.96754  .96625  .96507 

.96419 

.96369 

.10 

.96032 

.95726  .95396  .95088 

.94852 

.94714 

.15 

.94985 

.94421  .93796  .93200 

.92737 

.92461 

.20 

.93690 

.92771  .91734  .90727 

.89935 

.89459 

.25 

.92091 

.90698  .89104  .87538 

.86298 

.85550 

.30 

.90117 

.88106  .85785  .83493 

.81677 

.80582 

.35 

.87681 

.84886  .81647  .78454 

.75938 

.74433 

.  40 

.84679 

.80908  .76555  .72302 

.68995 

.67044 

.45 

.80980 

.76033  .70386  .64966 

.60843 

. 58460 

.50 

.76432 

.70115  .63052  .56463 

.51605 

.48878 

.55 

.70358 

.63025  .54541  .46946 

.41585 

. 38688 

.  60 

. 64068 

.54680  .44974  .36768 

.31299 

.28485 

.65 

.55885 

.45111  .34681  .26519 

.21472 

.19033 

.70 

.46204 

.34561  .24277  .17022 

.12953 

.11133 

.75 

.  35124 

.23630  .14697  .09215 

.06500 

.05398 

.80 

.23196 

.13419  .07076  .03845 

.02480 

.01981’ 

.85 

.11814 

.05482  .02310  .01049 

.00612 

.00474 

.90 

.03442 

.01150  .00357  .00128 

.00065 

.00049 

.95 

.00197 

.00037  .00006  .00001 

.00000 

.  00000 

50 


.000 

.99999 

.99998 

.99994 

.99985 

.99968 

.99935 

.99877 

.99780 

.99621 

.99370 

.98983 

.98405 

.97558 

.96347 

.94653 

.92338 

.89247 

.85216 

.80094 

.73767 

.66192 

.57438 

.47738 

.37515 

.27393 

.18134 
.10505 
.05049 
.01849 
. 00440 

.00045 
. 00000 


11 


\ 


TABLE  3 

CONFIDENCE  COEFFICIENT  ETA  WHEN  N*  50  AND  RT*.40 


WO  0.000  -0.500  -1.000  -1.500  -2.000  -2.500  -3.000 
RHO* 

.15  .99978  .99985  .99990  .99993  .99994  .99995  .99995 
.10  .99954  .99964  .99971  .99977  .99980  .99982  .99983 
.05  .99907  .99917  .99925  .99933  .99937  .99940  .99941 
.00  .99819  .99819  .99819  .99819  .99819  .99819  .99819 
.05  .99658  .99624  .99585  .99548  .99519  .99502  .99494 

.10  .99375  .99252  .99104  .98952  .98827  .98749  .98714 
.15  .98890  .98571  .98168  .97734  .97363  .97125  .97015 
.20  .98079  .97376  .96452  .95424  .94525  .93945  .93674 
.25  .96762  .95366  .93488  .91371  .89517  .88322  .87766 
.30  .94681  .92130  .88684  .84829  .81509  .79405  .78439 

.35  .91490  .87166  .81419  .75186  .70013  .66843  .65422 
.40  .86759  .79953  .71256  .62342  .55383  .51331  .49579 
.45  .80016  .70106  .58263  .47124  .39158  .34840  .33061 
.50  .70852  .57636  .43333  .31410  .23842  .20106  .18660 
.55  .59120  .43246  .28280  .17677  .11921  .09398  .08494 

.60  .45229  .28499  .15411  .07923  .04596  .03339  .02930 
.65  .30434  .15589  .06530  .02606  .01253  .00828  .00704 
.70  .16862  .06487  .01934  .00557  .00213  .00126  .00104 
.75  .06880  .01780  .00338  .00064  .00018  .00010  .00008 
.80  .01681  .00250  .00026  .00003  .00001  .00000  .00000 

.85  .00159  .00011  .00000  .00000  .00000  .00000  .00000 
.90  .00002  .00000  .00000  .00000  .00000  .00000  .00000 


TABLE  4 


CONFIDENCE  COEFFICIENT  ETA  WHEN  N*  50  AND  Rj,* 


W0  0.000  -0.500  -1 
RHO* 

.05  .99980  .99978 
.10  .99958  .99947 
.15  .99913  .99880 
.20  .99824  .99736 
.25  .99651  .99440 


000  -1.500 


2.000 


2.500 


99975 

99934 

99835 

99609 

99115 


99972 

99919 

99783 

99452 

98699 


99970 

99907 

99736 

99303 

98293 


99968 

99899 

99704 

99200 

98010 


.30 

.35 

.40 

.45 

.50 


.99322 

.98707 

.97580 

.95568 

.92079 


.98847 

.97704 

.95586 

.91838 

.85563 


.98092 

.96085 

.92387 

.86039 

.76035 


.97103 

.93966 

.88295 

.78991 

.65389 


.96130 

.91902 

.84440 

.72734 

.56744 


.95450 

.90481 

.81870 

.68773 

.51663 


.55 

.60 

.65 

.70 

.75 


.86259 

.77053 

.63547 

.45827 

.26263 


.75765 

.61803 

.44232 

.25695 

.10631 


.61852 

.44295 

.26144 

.11496 

.03200 


.48119 

.29766 

.14295 

.04737 

.00900 


.38279 

.20976 

.08562 

.02297 

.00334 


.33047 

.16873 

.06297 

.01516 

.00194 


.80 

.85 

.90 


.09975 

.01699 

.00047 


.02473 

.00200 

.00002 


.00425 

.00015 

.00000 


.00071 

.00001 

.00000 


.00019 
.00000 
.  00000 


.00010 
.  00000 
.00000 


.50 


.000 

.99968 

.99895 

.99689 

.99150 

.97871 

.95117 

.89794 

.80659 

.66974 

.49471 

.30936 

.15352 

.05542 

.01287 

.00160 

.00008 
. 00000 
. 00000 
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